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Abstract 

We present the results of a recently developed approach where the interplay between the 
itinerant and localized character of electrons in narrow band materials is described by adding 
on-site correlation effects to a first principle band calculation: the single particle band states 
are treated as mean field solutions of a multi-orbital Hubbard Hamiltonian and the many- 
body term associated to localized e-e interaction is described in a configuration-interaction 
scheme. The method allows to calculate hole and electron spectral functions which can be 
directly compared with spectroscopical results. 

Introduction 

The Hubbard model, dominated by the competition between inter-site hopping and on-site 
electron-electron repulsion, is believed to describe the physics of narrow band materials such 
as transition metals, transition metal oxides, cuprates, etc. ||^]. In these systems the itinerant 
character of valence electrons - clearly shown by the k-dispersion observed in photoemission - 
coexists with strong local correlations responsible of other spectroscopical features - satellites, 
band-narrowing, and opening, in some cases, of a Mott-Hubbard gap. 

In spite of the enormous amount of theoretical and experimental work which has been done on 
cuprates since the discovery of high Tc superconductors, an unified theoretical description of the 
whole valence spectrum, from the high binding energy region characterized by satellites, up to 
the valence band top, including both unperturbed single particle like and strongly correlated Cu 
derived structures, is still missing; this is due to the difficulty to combine an accurate treatment 
of many body terms with a realistic description of the band structure. 

Most of the work on this subject has been based on a drastic simplification of either the band 
structure or the e-e interaction; if the complex structure of these systems is reduced to a model 
description involving only CuO planes, the e-e interaction can be treated accurately, for instance 



by exact diagonalization techniques of finite (and small) CuO clusters |j|] , by refined mean- field 
solutions of a two dimensional t-J Hamiltonian ||2| etc. . However, in this search for the simplest 
model containing the relevant physics of superconductors one may miss some important effects, 
related for instance to the coupling between adjacent CuO planes and the possibility of a 
quantitative comparison with spectroscopical results. 

Photoemission data of highly correlated materials have been also interpreted using theoreti- 
cal approaches based on impurity and cluster configuration-interaction models which assume a 
strong wave function localization and adopt a rather simplified description of the band structure, 
with a considerable number of adjusting parameters. They have been widely used to describe 
the main structures and satellite peaks observed in the angle-integrated photoemission spectra 
of CuO and of cuprate superconductors . Other approaches have been proposed based on 
the density functional approximation (self-interaction corrected Q and LDA-I-U functionals) 
which fully include the itinerant character of electron states but describe the electron-electron 
interaction as a mean-field effective single-particle potential. 

A theoretical approach is then needed which includes both the hybridization between Cu and 
the ligands (or between sp and d states in the case of transition metals) accounted for by first 
principle band theory, and a treatment of e-e interaction which must be non-perturbative - to 
deal with systems which are in the high correlation regime - and beyond mean field - to include 
finite life-time excitations. 

The 3BS method ^, |lO|, ^] can be seen as an extension to the solid state of the configuration- 
interaction scheme used for finite systems: the Hubbard Hamiltonian is projected on a set of 
states obtained by adding a finite number of e-h pairs to the ground state of the single-particle 
Hamiltonian and this expansion is truncated to include one e-h pair. The effect of electron 
correlation on one electron removal energies from a partially filled band is then described as 
hole-hole and hole-electron interaction. The 3BS theory corresponds to the solution of a 3-body 
scattering problem involving two holes and one electron and has been originally formulated by 
Igarashi Self-energy corrections, spectral functions and quasi-particle band structure can be 
calculated for systems in different correlation regimes, getting a complete picture of the whole 
valence spectrum, including both long-lived coherent quasiparticle structures and incoherent 
short-lived ones. 



Hamiltonians 

Since we want to augment band theory with the inclusion of on-site correlation it is essential 
to define the relationship between band and Hubbard Hamiltonian in order to avoid double 
counting of e-e interaction. The exact many body Hamiltonian in second quantization is 
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Here eia and ti^j/s are the intra- and inter-atomic matrix elements of the one-particle Hamilto- 
nian and C/q^, Ja/3 are on-site Coulomb and exchange terms: 

Ua/3 — ^iacr,if3a,if3cr,iaa — ^iaa,if3—a,if3—a,ia(Tj 

with 

ss' ' ' 

Different approximations to the exact Hamiltonian ([l|) can be obtained using a mean field 
approach which amounts to neglect fluctuations in the electron occupation 

+ {fliacj - (fliaa)) {flifScT' " {nif3a')) 

The mean field approximation can be applied to all the many body terms of (||) transforming it 
into a single particle Hamiltonian 

jjMF ^ ef^^hiaa + X! X! (Sclaa^j f^a , (2) 

or selectively to the multi-center integrals, keeping the full many body character in the one-center 
terms; in this way one gets the generalized Hubbard model 
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Since H^^^ and differ only for the treatment of the on-site correlation - included in 
as a mean field and treated as a many body term in - it is easy to show that 

^Za = €a + E [(f^"/3 - + ^afS {nif3~a)] ■ (4) 

Using a Bloch basis set the two approximate Hamiltonians become 
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here Oko-'^kl- ^'■^ destruction/creation operators of electrons with wave vector k , spin o", band 
index n and C"(kiT) are the expansion coefficients of Bloch states in terms of localized orbitals. 
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The relationship between single particle eigenvalues and e? appearing in the two Hamil- 
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which is the analogue of eq. (^) for Bloch states. Any band structure calculation corresponds 
to the solution of some H^'^^ describing the interacting system as an effective single particle 
problem and equations contain the correct recipe to include Hubbard correlation starting 
from band structure eigenvalues, avoiding e-e interaction double counting. 

Hole spectral function, self-energy and the Faddeev method 

We are interested in the hole spectral function 

n 

which describes the removal of one electron of wave-vector k, band index n and spin a and is 
related to the hole-propagator 

G-{\^na,uj) = -(^olal^lG(^)aLI'fo), z = -uj + E^{N,) +\6- 

(10) 

£^0 (-^e) and j^'o) define the ground state of the A'^e particle system and 

z - 

is the resolvent operator. By projecting the Hamiltonian (||) over a complete set appropriate for 
the Nf> — \ particle system one gets an expression for appropriate to describe one electron 
removal. The key approximation is to choose a subset of all the excited states of the non- 
interacting system and assume it to be complete. Any N particle non-interacting state can be 
obtained by repeated applications of creation/destruction operators to the ground state Slater 
determinant |<l>o) i.e. by adding e-h pairs to it; according to 3BS the interacting state with one 
removed electron of momentum k and spin a is expanded in terms of the basis set including 
1-hole and 3-particle configurations 

\s) = akn<T l^o) , I*) = «q3n3cr3«q2n2fT2«qiniCTi \^o) , (12) 

with 

qi + q2 - qs = k, ai 0-2 - 0-3 = a. 

The effective Hamiltonian for the N-1 particle system is then 

^H^+H^ + V, (13) 

where H\ is associated to one-hole configurations 

= {s\H^\s)\s) 
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describes the contribution of 3-particle configurations 

H, = Y^{t\H"\t')\t){t'\, 
w 

and V is the couphng between 1- and 3-particle states 

V = Y,{s\H^ \t)\s) {t\ + h.c. 



We can now calculate the resolvent (11). We define the 3-particle resolvent, that is the resolvent 
associated to the 3-particle interaction 



and the Dyson equation which relates G{z) to it 

G{z) = h{z) + h{z)[H^ + V\G{z). (14) 

The Faddeev scattering theory allows to determine F2,{z) by separating the 3-body Hamiltonian 
in diagonal and non-diagonal parts 



H^ = Y.{t\H''\t)\t) {t\, 
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Hr=T^{t\H''\t') \t){t'\ 
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defining the diagonal 3-body resolvent 



S = Hi:"' + H?"fPS^ 



and the scattering operator 

^-^'3 ^-3 -3 

The full 3-body resolvent can be written in terms of the diagonal one and of the scattering 
operator as 

F, = Fi> + Fi>SFi>. (15) 
As shown in references ||9|, [l^ the non-diagonal 3-body interaction is the sum of two potentials 

H^^ = Vh-h + Vh-e, 

which describe h-h and h-e multiple scattering. We define partial scattering operators Sh-h,Sh-e 
such that S = Sh-h + Sh^e-, i-e. 

Sh-h = Vh-h + Vh-hFf S , 
Sh~e = ^h-e + ^h-eFf S , 

which are related to the scattering T-matrices 

fh-h = Vh-h + Vh-hFi'n-h, (16) 

fh-e = Vh-e + Vh-eFi'th-e, (17) 




Figure 1: Pictorial representation of non-diagonal terms in the effective hole Hamiltonian: V 
couples one- and three-particle configurations while Vh-h and Vh-e describe scattering between 
three-particle states, namely hole-hole and hole-electron scattering respectively. 



through the Faddeev equations |12| 



Sh-h = Th-h + Th^hFf Sh-e, 

Sh-e = Th-e + Th-eF^Sh-h- (18) 

Inserting (18) into ( p^ ) one gets the expression for the 3-particle resolvent in terms of scattering 
operators S^-e and T^-^ 

#3 = + Ff {fh-h + fh^hFfSh^e + Sh-e) F^ . (19) 

In this expression and T-operators -or rather their matrix elements between three-particle 
states- have an analytical expression, while the inclusion of scattering operator S^-e will require 
the solution of an integral equation. 
After some algebra the hole propagator becomes 

g- (kna, uj) = -Gss {z) = \ — ; (20) 

w 

with the notation Gss = {s\ G |s), F^tt' = {i\ F3 \t') etc.. Since the difference between the ground 
state energy of the A'^g-particle system Eq{N(,) and the average of over \s) states turns out 
to be 

Eo{Ne) - = e^^^ + = e^f^, 
the mean field band eigenvalues appear naturally in the denominator of the hole propagator. 



Comparing eq. (20) with the usual expression 
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we can identify the self-energy correction to band eigenvalues as 

E- (knc7, u) = -J2 Fm' VfsVst. (22) 

tt' 

The procedure we have outlined ends up with a result which has a simple physical interpretation: 
the creation of one hole in an unfilled valence band is followed by multiple h-h and h-e scattering 
which is responsible of a renormalization - through self-energy corrections - of the energy states. 
The efficiency of the scattering processes depends a) on the strength of the screened on-site e-e 
interaction and b) on the available empty states (i.e. on the number of initial valence holes). 
This explains the well known differences between various transition metals (Ni and Cu, for 
instance) and possibly those arising in cuprates as a consequence of hole doping. 
The self-energy S"^ and spectral function for electron addition can be calculated in the same 
way as described above just exchanging empty states with filled ones and vice versa, as described 
in detail in ref. @. 

3BS at work 



In order to calculate the self-energy S~ (kna,uj) according to (22) one has to perform summa- 
tions over the 3-particle states \t) involving k-vector conservation; this is done within the so 
called local approximation 

ix 

This approximation allows to transform the k-vector summations into integrals involving the 
orbital density of states nj^e) and to calculate T-matrices, scattering operators, resolvents, self- 
energy and finally spectral functions according to the computational strategy discussed in detail 
in ref. |9|, ^ which can be summarized as follows: 

• Input: band structure (e^, C2|(fe), na|(e)) and U; 

• free propagators 



9h-hi^> 



oo J —oo ^ 



e' — e — i6^ 



J_oo JEf uj-e' + e-io J-oo uj-e'-io 



T-matrices 



kernel 

6, e') = [\e"n^,{e")g^u; + e" - e)g\u: + e" - e')Tf_^{u: + e")Tf_^{u: - e"), 



and 

rE 



B^\u,e) = /■ deVi(e')/(^ + e'-e)r,"f^(u; + e') 

^— oo 

roo 

gf.S^ -e')+ de"n„i(6")/(^ + e' - e")gl\{u - e")Tf_^{u^ - e") 

JEf 
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Figure 2: Real (solid line) and imaginary (dotted line) part of self-energy T,^^{uj) (a) and spectral 
function (b) in paramagnetic CuGeOs for the creation of a hole at the T point corresponding to 
the single particle band energy = -3.58 eV. Energies are referred to Ej. 



solve the integral equation 



POO 

A"'3(a;,e) = 5°^(a;,e) + / d eVi(e')^"^('^, e'); 

JEf 

orbital self -energy 

^^t(^) = E L den,i(6)r^f,(u; - e) [l + UA-^{u;e) 



K- and band-index dependent self-energy 



^~{kn],u;) = Uj2\CJ^^{k)\' 





/ den„|(e)S^^(u;) 

„, J—oo 



spectral function 



d\ \u;) = ^ Im 

^ kn 



uj — — Tj {kn |, w) 



We report the results for a transition metal (nickel) and a cuprate (CuGeOa) obtained by 
considering the interaction between opposite spin electrons localized on the transition metal 
sites as the dominant contribution, i.e. 

Udd for a, /3 = d orbitals 
. elsewhere; 

Ua(3 - JcP ^ 0. 



Ual3 



To apply this method to CuGeOs we have used the eigenstates/eigenvalues of ref. [14] and 
assumed Udd = 8 eV. Fig. |2| shows the hole spectral function and self-energy for a particular 
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Figure 3: Total spectral function D(u;) = D {to) + D+(a;) for CuGeOs calculated by (a) 3BS 
theory with U=8 eV; (b) single particle density of states (Mattheiss, PRB 49, 14050 (1994)). 
Ej corresponds to lij = eV. 



eigenstate (ej^=-3.58 eV at the F point). The peaks in the spectral function can be classified 
either as quasiparticle excitations or as satellites, according to the value of the imaginary part of 
the self-energy in the region of the peak: quasi-particle excitations correspond to small imaginary 
part, and give rise to the coherent part of the spectral function. Satellites occur where the 
imaginary part of self-energy is large and correspond to short-lived excitations with a large 
intrinsic line-width; we refer to them as to the incoherent part of the spectral function. 
The hole and electron spectral functions are plotted in fig. ^ and compared with single particle 
results. The effect of electron correlation on single-particle states is dramatic: some bands 
are shifted to higher binding energies, spectral weight is removed from the upper part of the 
spectrum, and many new states (satellites) appear; only states around —8 -. — 9 eV and —5 -i- 
—6 eV are practically unaffected being mainly Ge and O derived. CuGeOs is an insulator but it is 
predicted to be a metal by single particle band calculation; the inclusion of electron correlation 
reproduces this insulating behavior and the energy gap, calculated as the energy separation 



between electron removal and electron addition spectra, reproduces the experimental one [15|. 
The same was proven to be true also in the case of NiO, where a 3BS description of Hubbard 
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Figure 4: (a) Energy position of the peaks of nickel in the spectral function at k-points along 
high symmetry lines of the Brillouin Zone, (b) Single particle band structure (Manghi et al., 
PRB 56, 7149 (1997)). 

correlation was able to reproduced both the complex satellite structure and the measured value 
of the insulating gap [p^]. 

The ability of 3BS approach to open up Hubbard gaps, i.e. to reproduce an insulating behaviour 
in a system which is metallic according to band theory, is related to its non-perturbative char- 
acter; the method can be applied in any correlation regime and it has been shown [^, 17| 



to reproduce for U much larger than the band width W {U^ oo) the so called "Hubbard I" 



solution |18| of Hubbard Hamiltonian, i.e. the atomic limit where hole and electron states are 
separated by a Mott-Hubbard gap equal to U [|l9|] . 

We describe now the results of the application of 3BS to nickel. In this case the on-site e-e 



repulsion is more effectively screened and the estimated value of Udd is — 2 eV |20]. Fig. 



reporting the comparison between quasi-particle states and single particle ones, shows that e-e 
correlation effects are still sizable and they are actually essential in order to to reproduce the 
observed spectroscopical features, i.e. satellite structure at 6 eV binding energy, correct band 
width (overestimated in LDA), exchange splitting |2^, and energy dispersion |1C]. 



Summary and outlook 

We have described a method to include on-site interactions in the description of hole and elec- 
trons states: ab-initio single-particle band states are used as input mean-field eigenstates for the 
calculation of self-energy corrections according to a 3-body scattering (3BS) solution of a multi- 
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orbital Hubbard Hamiltonian. When applied to valence states of ferromagnetic nickel it allows 
to get a quasiparticle band structure which compares much more favorably with the experimen- 
tal observation than conventional mean-field LDA, reproducing the observed band width, the 
energy dispersion, the satellite structure and the exchange splitting. Since the method does not 
rely on a perturbation expansion it has a wide range of application, including any correlation 
regime. In the case of a highly correlated system such as CuGeOa, 3BS is able to reproduce 
both the insulating behaviour and a correct overall picture of photoemission experiments. 
Our present choice of empirically determining the parameter U of the Hubbard Hamiltonian - 
which has been fixed to reproduce the satellite binding energy - ensures that we obtain a good 
agreement with experiments; however the possibility of reproducing both the satellite structures 
and other spectroscopical features such as energy dispersion and spin dependence in nickel and 
the insulating gap in CuGeOs can be seen as a non trivial result and a success of the method 



itself: previous methods based on a simplified description of the scattering channels p2| , 23 1 
in fact have not been able to reproduce at the same time the satellite energy position and 
the valence band width which turned out to be systematically overestimated for values of the 
Coulomb integral fixed to reproduce the satellite binding energy. The problem of extracting 
Hubbard U from ab-initio calculation, either in the so called Constrained-Density Functional 



scheme [ p4[] or as screened Coulomb interaction 25|, is an important issue which goes in 
the direction of a full match between model Hamiltonians and realistic systems and that we are 
presently considering as an implementation of our approach. 
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